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Abstract. It is known from both experiments and molecular dynamics simulations 
that chemically patterning a solid surface has an effect on the flow of an adjacent 
liquid. This fact is in stark contrast with predictions of classical fluid mechanics 
where the no-slip boundary condition is insensitive to the chemistry of the solid 
substrate. It has been shown that the influence on the flow caused by a steep 
change in the wettability of the solid substrate can be described in the framework 
of continuum mechanics using the interface formation theory. The present work 
extends this study to the case of intermittent patterning. Results show that vari- 
ations in wettability of the substrate can significantly affect the flow, especially 
of thin films, which may have applications to the design of microfluidic devices. 

1 Introduction 

Flows of liquids on small length scales, characterised by high surface-to-volume ratios, are of 
increasing interest in many emerging technologies [1 . The correct description of the physics 
at liquid-solid interfaces is imperative to successful modelling of such phenomena [2]. Partic- 
ular interest lies in how modification of the chemical properties of the solid substrate, most 
importantly variation in its wettability, may affect the flow of an adjacent fluid. In classical 
fluid mechanics, the no-slip boundary condition is applied on a solid substrate irrespective of 
its chemical properties, whereas molecular dynamics simulations [3] and experiments [4] have 
shown that deviations from this condition can be significant, even for chemically homogenous 
substrates, when considering flows on small length scales. The conventional generalisation of 
no-slip is the Navier slip condition [5] , which states that slip (a non-zero difference between the 
tangential velocities of the fluid and substrate) is proportional to the shear stress acting on the 
liquid-solid interface from the fluid. Using the coefficient of proportionality as an adjustable pa- 
rameter, one can describe flow over chemically homogenous solid substrates in situations where 
slip is important (e.g., [6]). However, the Navier condition is a dynamic boundary condition and 
as such it is unrelated to the solid's wettability, which can be defined in the absence of motion, 
for example in terms of the static contact angle between a free surface and the solid substrate. 
Consequently, the Navier condition offers no framework that would allow one to incorporate 
the solid's wettability and, in this sense, has no advantages over the no-slip condition. 

It has been shown [7] that the interface formation theory [8] that has previously been 
applied to describe dynamic wetting and some other phenomena |9ll0lllj can be used, without 
any alterations, to study the effect of a variation in wettability of the solid on the flow. In 
the present work, we consider a particular application of the theory in the situation where the 
pattern on the solid substrate consists of an array of stripes of differing wettabilities. 
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Fig. 1. Problem setup. 



2 Problem formulation 

Consider the steady two-dimensional flow of an incompressible Newtonian liquid between two 
solid surfaces (Fig. [1]) where the upper surface is chemically homogenous and moves parallel 
to itself with velocity u = (u, v) = (U, 0), whilst the lower surface is chemically patterned and 
remains stationary. The pattern on the lower surface consists of intermittent stripes of solid 
1, with width L — a, and solid 2, with width a. It is assumed that both surfaces are perfectly 
smooth. The wettabilities of solid 1 and solid 2 are characterised by the contact angles, 0i 
(i = 1,2), which a liquid-gas free surface would form with solid i. 

In the bulk the flow parameters satisfy the Navier-Stokes equations, 

V-u = 0, pu • Vu = -Vp + /iV 2 u, (1) 

where p is the pressure, p is the density and p is the viscosity of the fluid. On both surfaces we 
apply the boundary conditions provided by the interface formation theory [8]: 

<* = l(P(0) - P% (2) 
o T + TT )=0u, pv = — , (3) 
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v =* u + a \Tx + ^^)- (5) 

A detailed explanation of the process of interface formation together with a derivation of the 
equations and applications of the theory are given in several works (e.g. [819110] ) so that here 
we will briefly recapitulate only the important points. 

At constant temperature the state of the liquid-solid interface, i.e. physically, a microscopic 
layer of liquid adjacent to the solid surface, is characterised by the surface density p s of liquid 
in this layer. In equilibrium one has p s — p s e and the surface tension a{p s e ) — cr e , where the 
equilibrium values p s e and a e are different for solids of different wettability. The simplest equation 
of state accounting for the influence of p s on a is given by ([2]). Instead of no-slip one now has 
(J3j) , where the first equation is the generalised Navier condition given here for the case of a 
variable pi and the second describes the mass exchange between the bulk and the liquid-solid 
interface. This mass exchange is also accounted for in the mass balance equation for the interface 
(j4|). Finally, equation (O describes how the velocity in the interface v s depends on the bulk 
velocity u evaluated on the liquid-facing side of the liquid-solid interface and the surface tension 
gradient. The parameters a and (3 characterise the response of the interface to surface tension 
gradients and external torque, respectively; 7 is associated with the inverse compressibility of 
the fluid; p* ^ is the surface density corresponding to zero surface tension; t is the surface 
tension relaxation time; in the simplest variant of the theory these are all treated as material 
constants. 

The terms including dp s e /dx in ([3]) and ([5]) describe the surface force which ensures that, 
in the case of a variable p s e (and hence cr e ), there is no perpetual motion. Once the outer flow 
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drives the matter in the liquid-solid interface along a solid surface whose wettability varies, the 
surface tension gradient in ([3]) and ([5|) becomes unbalanced triggering the interface formation 
process which in its turn affects the flow that caused it. 

The last step now is to specify the wettability pattern of the lower surface, i.e. prescribe 
the dependence of p 3 e on coordinates along it and the value of p 3 , on the upper surface . From 
the viewpoint of fluid mechanics, the junction between solid 1 and solid 2 is associated with 
a discontinuity in the properties of the substrate and, therefore, requires special handling. In 
order to avoid this, we consider a transition region in which the equilibrium surface density 
varies smoothly, from its value p 3 e on solid 1 to value p| e on solid 2 (or vice- versa), across a 
finite distance This can be modelled, for example, using 

Pi = P\e + I (PL - Pie) {tanh[(a/2 + x)/l] + tanh[(a/2 - *)/*]} . (6) 

On the upper solid we assume p s e = p\ e . 

Using the Young equation, cr(p| e ) = —<Ji g cos where o\ g is the surface tension of the liquid- 
gas free surface, and the equation of state ([2]) one can relate p\ e (or p| e ) with the equilibrium 
contact angle Q\ (or 6*2) which the free surface would form with solid 1 (or 2): 

Pie = P(o) + 7 _1 ^ 9 COS0i, (i = 1, 2). (7) 

The contact angle is the conventional and convenient measure of the solid surface wettability 
with respect to a given fluid, and in what follows we use Bi instead of pf e . 

Due to the periodicity of the problem it is convenient to use periodic boundary conditions 

du . du 

u U=-L/2= u U=L/2> -q^ \x = -L/2— -q^ U=i/2) (°) 
P S |x=-£/2= P S \x=L/2, V s \ x= - L /2= V s \ x=L/2 ■ (9) 

Equations ([l])-©, dSJ) — d9j) now fully specify our problem. 



3 Results 

The problem has been solved numerically using the finite element method. After using L, 
d, U, dU/L, pUL/d 2 , ui g and p 3 ^ as scales for the horizontal coordinate, vertical coordinate, 
tangential velocities, vertical velocity, pressure, surface tension and surface density, respectively, 
one has that the problem is specified by following non-dimensional parameters: 

d _ a „ pUd uU Ut - (3Ud 
r = —, a = —, Re = , Ca = , e = — , f) = , 

L L P (Jig d (Tig 

Q = %a = ^,T=l,X= 7 -?k,e i (i = l,2). 
pUT Ud L aig 

As the interface formation theory is applied here without any ad hoc alterations, we may take 
estimates for the theory's phenomenological constants from independent experiments. Dynamic 
wetting experiments [llj . for example, suggest that for a fluid with p, ~ 10 g cm -1 s _1 and 
p ~ 1 g cm~ 3 we have a ~ 10~ 8 g _1 cm 2 s, p ~ 10 s g cm -2 s _1 , p 3 ^ ~ g cm~ 2 and 

t ~ 10~ 7 s. Using these estimates and taking typical values associated with the flows of thin 
liquid films, U ~ 10 cm s _1 , d ~ 10 -5 cm, ai g ~ 10 2 dyn cm -1 and !~5x 10~ 7 cm, we have 
the following typical magnitudes of the non-dimensional groups: Re = 10~ 5 , Ca = 1, e = 10 _1 , 
P = 10 2 , Q = 1Q- 1 , a = W 2 , 1 = 5 x 10" 2 and A = 2. 

Consider each stripe of solid to have equal width, a = 1/2, and the aspect ratio r = 1. 
Then, without loss of generality we may consider solid 1 to be more hydrophilic than solid 
2. Streamlines for the typical values and 9\ — 10°, 62 = 80° are shown in Fig. [2l The less 
hydrophilic substrate appears to acts as an obstacle to motion with streamlines coming out 
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Fig. 2. Streamlines and the corresponding components of velocity on the substrate for flow over a inter- 
mittently chemically patterned surface in which values of the streamfunction ip are given. Parameters 
are at their typical values with 81 = 10° and 62 = 80°. 



of the surface as the solid becomes less hydrophilic, at x = —0.25, and entering again as 
the solid becomes more hydrophilic again, at x — 0.25. This effect can be seen clearly by 
looking at the normal velocity on the liquid facing side of the liquid-solid interface in Fig. [5] 
One can see that the components of velocity are trying to relax to their equilibrium values, 
u e = (9.8 x 10~ 3 , 0) but are unable to obtain them before the next transition in wettability. 
This means that the distance between transitions is of the same order as the relaxation length, 
which is the distance that it takes for surface properties to relax to their equilibrium values. 
The width of the transition region I is far smaller than the relaxation length and, as in [7], its 
exact size has very little effect on the overall flow field. 

Notably, in contrast to studies of slip on super hydrophobic substrates, where the presence 
of nanobubbles leads to the effective slip length of the system becoming a relevant measure of 
the effects of chemically patterning the substrate [12|13j , here the slip coefficient (3 remains the 
same on both solids. Importantly, as shown in [7j, slip results primarily from the disturbance of 
the force balance in the liquid-solid interface and not from the tangential stress, as follows from 
the standard Navier condition. This has significant consequences when attempting to interpret 
the results of molecular dynamics simulations and experiments in terms of a continuum theory. 
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